Exactly solvable Kitaev model in three dimensions 
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We introduce a spin-^ model in three dimensions which is a generalization of the well-known 
Kitaev model on a honeycomb lattice. Following Kitaev, we solve the model exactly by mapping it to 
a theory of non-interacting fermions in the background of a static Z2 gauge field. The phase diagram 
consists of a gapped phase and a gapless one, similar to the two-dimensional case. Interestingly, 
unlike in the two-dimensional model, in the gapless phase the gap vanishes on a contour in the 
k space. Furthermore, we show that the flux excitations of the gauge field, due to some local 
constraints, form loop like structures; such loops exist on a lattice formed by the plaquettes in 
the original lattice and is topologically equivalent to the pyrochlore lattice. Finally, we derive a 
low-energy effective Hamiltonian that can be used to study the properties of the excitations in the 
gapped phase. 

PACS numbers: 75.10. Jm,03.67.Pp,71.10.Pm 
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I. INTRODUCTION 

The study of topological phases has been actively pur- 
sued in condensed-matter systems for some years. This 
has resulted in the emergence of a new paradigm in the 
theory of quantum phase transitions: certain phase tran- 
sitions cannot be described in terms of local order pa- 
rameters associated with spontaneously broken symme- 
tries; instead, the phases in such transitions arc char- 
acterized by topological order i The most famous ex- 
ample of topological order in a quantum system is in 
the phenomenon of fractional quantum Hall effect i^i^i^ 
Other examples — experimental and theoretical — include 
quantum spin liquids (^'^'^i'^ quantum dimer modelsf^i^iii 
quantum loop modelsfi^ etc. 

Recently, it has been proposed that topological phases 
can be used to do quantum computationj ^'^i^^ The main 
obstacle in the realization of quantum memory — the ba- 
sic ingredient of a quantum computer — is dccoherencc: 
it is difficult to prepare states that arc robust to exter- 
nal noise. Kitae^ii^ suggested that topologically ordered 
states can be used to overcome this problem. He illus- 
trated these ideas in a spin-^ model on a hexagonal lat- 
tice which, quite remarkably, can be solved exactly. The 
ground state of the Kitaev model has two phases: in one 
phase the elementary excitations have a gap in the spec- 
trum and are Abelian anyons; the second phase is gapless 
in the absence of an external magnetic field, but develops 
a gap when the field is switched on, and then the exci- 
tations are non-Abelian anyons. In topological quantum 
computation, braiding of non-Abelian anyons is essential 
for the realization of universal quantum gates Regard- 
ing the feasibility of a physical realization of the Kitaev 
model, there has been a proposal to realize it on an op- 
tical latticeJiii^ 

Apart from its potential application in quantum com- 
putation, the Kitaev model is an interesting many-body 
system by itself. First, exact solutions are rare in dimen- 
sions higher than one, and second, the model provides 



a relatively simple platform — the Hamiltonian involves 
only two-body interactions — to study concepts such as 
topological order and fractional excitations. Thus, 
not surprisingly, the various many-body aspects of the 
model have been thoroughly investigated . ^^^^'-''^^i^^'^'^i^^ 
There have also been some generalizations to other two- 
dimensional (2D) latticeS i^^'^^ It is then worthwhile to 
find generalizations of the model in higher dimensions. In 
this paper, we introduce and study a three-dimensional 
(3D) version of the Kitaev model. Three-dimensional 
models exhibiting topological order have been studied 
previously.i22i2^'^'^ 



An exact solution of the Kitaev model is possible due 
to the existence of a macroscopic number of locally con- 
served quantities: this facilitates the mapping of the 
model to a quadratic Hamiltonian of Majorana fermions 
hopping in the background of a static Z2 gauge field. 
The three-dimensional model we construct also has the 
above feature, which renders it exactly solvable. It has 
a gapped phase and a gapless one, just as in 2D, and 
although the phase boundaries are identical to the lat- 
ter, the nature of certain excitations is quite different. 
In the three-dimensional model, the excitations of the 
gauge field are localized on "loops" . This gives rise to 
the possibility that, as yet unverified in our model, such 
excitations can obey nontrivial statistics, since double 
exchange in loops — unlike points — is topologically non- 
trivial in three dimensions (for example, see Ref. \2w . 

The paper is organized as follows: in Sec. |TT]we define 
the spin Hamiltonian and then rewrite it using a Majo- 
rana fermion representation of spin i. In Sec. IIIII the 
low-energy spectrum for fcrmionic excitations is derived 
while Sec. IIVI discusses the excitations of the gauge field. 
We then derive the low-energy effective Hamiltonian in 
the large Jz limit in Sec. FVland end with a discussion of 
our results in the last section. 
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FIG. 1: The honeycomb lattice: the three types of hnks are 
labeled x, y, and z. 



II. 



HAMILTONIAN 



Before presenting the three dimensional model, let us 
briefly describe the Kitaev model in 2D. It is a system 
of spin-i degrees of freedom located at the vertices of a 
honeycomb lattice. In the honeycomb lattice, there are 
three types of links which, distinguished by their orien- 
tations, are labeled x, y and z (see Fig. [1]). The Kitaev 
Hamiltonian is 



y E 

<l,3>y 



J. 



E 



<i-3>z 



(1) 

where cr^'s are the Pauli matrices, and < i, j >a indicates 
that i and j belong to a link of a-type. 

For our purposes, the key points to note about the 
Hamiltonian in Eq. ^ are that for any spin, only one 
of the components couples to a particular neighboring 
spin, and the component with non-zero coupling strength 
is different for each of the three neighbors. As we will 
see later, this leads to the existence of a set of mutually 
commuting conserved plaquette operators, which in turn 
makes the problem exactly solvable. 

We note two features of the honeycomb lattice that are 
pertinent to the construction of the 3D lattice: (1) The 
coordination number of the lattice is 3. (2) The three 
types of links x, y and z are distributed in such a way 
that two links of the same type do not touch each other. 



A. Lattice 

To facilitate visualization, we will first describe how 
to obtain the 3D lattice starting from the familiar cubic 
lattice. Let i,j,k S Z be the x,y and z coordinates 
of the latter. The new lattice is obtained by removing 
those sites that satisfy one of the following conditions: 

(1) fc = mod 4 and i = mod 2, (2) fc = 1 mod 4 and 
j = mod 2, (3) k = 2 mod 4 and i = 1 mod 2, and (4) 
k = 3 mod 4 and j = I mod 2. 

This amounts to depleting the cubic lattice by half, and 
the resultant lattice has coordination number 3 (see Fig. 

(2) . We note that: (i) The x-y planes alternately consist 



of disconnected rows or disconnected columns, (ii) As 
one goes along a particular row (column), at each site 
there is a link whose direction alternates between positive 
and negative z axes. That is, there is a criss-crossing 
structure between adjacent planes which ensures that the 
lattice is truly three dimensional — despite a coordination 
number of 3 — and not a set of mutually disconnected two- 
dimensional surfaces. 
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FIG. 2: The 3D lattice: the four sites inside the loop (marked 
1-4) constitute a unit cell; ai , a2 , and as are the basis vectors. 
Plaquette p consists of sites marked 1-10. 



To parametrize the lattice sites, we first note that the 
unit cell contains four sites. The position vector of a unit 
cell is given by 

R = msii + na.2 + pEL3 , m,n,p'El,, (2) 
ai = 2x, a2 = 2y, a,, = x + y + 2z, (3) 

where x, y, and z arc unit vectors along x, y and z direc- 
tions, respectively. The four sites within a unit cell are 
at 



ri = R 



y 

2 



z, r2 



R 



y 

2' 



r3=R+-, r4 = R+- 
^ 2' 2 



(4) 



To define a Kitaev-type Hamiltonian, we need one more 
ingredient, viz, the labeling of links. To this end, we al- 
ternately assign x and y labels to the links in each of the 
rows and columns that lie on the x-y plane; the remain- 
ing links, the ones along the z axis, are labeled z. [The 
ambiguity in the assignment of x and y labels within each 
row (column) is resolved by demanding periodicity.] This 
way of labeling ensures that the three links emanating 
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from each site have different labels. Now the definition 
of the Hamiltonian in Eq. ^ can be applied to the three- 



dimensional lattice we have constructed. Explicitly, 



J 



H = ^ [- J,a^(R)af(R-a3)- Vf(R)af(R + ai-a3)- J,af(R)af(R) 

R 

- Jyaf (R)(Tf (R + as) - J^cr^ (R)ct|(R) - J,al{R)al{Ii) 



(5) 



r 



For simplicity of notation, we will continue to use the for- 
mal expression for H in Eq. ([T]), where sites are referred 
to by a single lower-case index, and revert to the explicit 
form in Eq. ([5]) only when the calculation demands it. 



B. Conserved quantities 

Let I = (ii7 j2, • ■ • , jn) be a sequence of lattice sites 
such that jm and jm+i are neighbors for to = 1, 2, ... n 
with jn+i = ji- Then I represents a loop. Equivalently, 
the links formed among jm also uniquely determine the 
loop. Let 



W, 



(6) 



where is the label of the link that connects jm to its 
neighbor which is not in /. It is easy to see that 

[Wi,H]=0, [Wi,Wi,]=0, yij'. (7) 

However, not all such operators are independent. To 
see this, we define a geometrical operation of combining 
two loops as follows: the combination of loops and I2 
is the loop li2 formed by all the links in li and I2 ex- 
cept those that are common to both. Then, it is easy 
to see that Wi^Wi^ = -izWi^^. This means that for a 
lattice with open boundary conditions, any loop opera- 
tor can be written as a product of those defined on the 
plaquettes — the elementary loops which cannot be ob- 
tained by combining smaller loops. For example, in Fig. 
m Wp corresponding to the plaquette p, consisting of the 
sites marked 1,2,..., 10, is 



(8) 



There are further constraints among Wp corresponding 
to different plaquettes, and we will discuss them in Sec. 
TVl where the excitations are studied. 



Hamiltonian in Majorana fermion 
representation 



To solve the Hamiltonian, following Kitaevfi^ we use 
a representation of Pauli matrices in terms of Majorana 



fermions. At each site j, we introduce four Majorana 
operators 5J, 6^, 6| and Cj, which satisfy the following 
relations. 



S ~ ^3 ' 



(9) 



The Majorana fermions act on a Fock space which is 
four dimensional, whereas the spin Hilbert space has two 
dimensions. We define the physical space as consisting 
of those states which satisfy the constraint 



^^■10 = 10, Vj where i?, = 6|6^^6|r 



(10) 



= 1; thus, {l,Dj) form the elements of a Z2 gauge 
group. The spin operators are defined as 

= ih^^c,. (11) 

When restricted to the physical space, the above opera- 
tors satisfy the standard spin-i algebra. The projector 
to the physical space is given by 



(12) 



All the states related by a gauge transformation project 
on to the same physical state. 

Starting from a generic spin model, i/jcr"}, we can 

obtain a fermionic Hamiltonian, Cj}, by using Eq. 

(fTTj) . Since \H,Pj\ = Vj, the eigenstates of H can be 
obtained from those of H by projecting the latter to the 
physical space. 

Substituting Eq. ((ll|) in Eq. H]), we obtain 



A 



jk 



Ujk 



Jotjk^jk if j and k are linked. 







— '"i "k ' 



otherwise, 



(13) 

(14) 
(15) 



where ajk is the type of the link between j and k. We 
note that Ujk = —Ukj, and in the sum the links are 
treated as directed and therefore counted twice. We use 
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a hat to emphasize that ujk is an operator; Ujk is the 
corresponding eigenvalue and takes vahies ±1. Further- 
more, 

[H,Ujk]=0, [ujk,uir,i] = Oy < j,k >,< l,m > . (16) 

Therefore, the Hilbert space breaks up into various sec- 
tors, each corresponding to a particular set {ujk}] the 
matrix elements of H between states belonging to differ- 
ent sectors are zero. The Hamiltonian in a given sector 
is obtained by replacing the operators Ujk with their cor- 
responding eigenvalues, Ujk'- 

Hu^^^AjkCjCk, (17) 

where Ajk is obtained from Eq. (|14[) by substituting 
Ujk with Ujk- The gauge-invariant and hence physical 
conserved quantities are 

_ 10 
m— 1 

where p = {ji,j2, ■ ■ ■ ,jio); as before, is a plaqucttc. Wi 
is related to Wi in Eq. ^ as follows: 

Wi = riWiri, ri = l[Pj. (i9) 

From now on we will simplify the notation, following 
Kitaev^ by not making distinction between operators 
in the physical and extended Hilbert spaces, i.e., we will 
drop tilde from all the operators acting in the latter. 

III. GROUND STATE AND SPECTRUM 

Next question to be addressed is: Which sector of 
{uij} does the ground state belong to? The problem 



of free fcrmions hopping on a d-dimensional hypercubic 
lattice with hopping amplitude |iy |e*^'J between nearest- 
neighbor sites i and j was studied by Lieb.— Moreover 
it was shown that, if |ty| is reflection symmetric about 
certain planes that does not contain any sites, then the 
ground-state energy is the lowest when the flux of the 
phase along the plaqucttcs, <!> = X]<y> ' equals tt if 
the length of the loop is mod 4, and zero if the length is 
2 mod 4. Unfortunately, Lieb's result cannot be directly 
applied to our case because our lattice does not have the 
required reflection symmetry. 



We numerically studied the ground-state energy of the 
Hamiltonian in Eq. pT)) for lattices containing up to 864 
sites using periodic boundary conditions. The ground- 
state energy in the sector with uniform flux (Wp = 
+l,Vp) was compared with that of other selectively cho- 
sen flux configurations. As will be discussed in Sec. llVi 
due to some local constraints, the plaqucttcs that are ex- 
cited themselves form loops in an embedded lattice. The 
smallest such loop consists of six plaquettes. We consid- 
ered excitations of loops of varying length and found the 
energy increasing with increasing length. We also looked 
at excitations of multiple loops. In every case considered, 
we found that the energy of the lowest energy state is 
greater than that of the Wp = +1 uniform flux state. In 
Sec. |Vl wc will get further confirmation of this result, at 
least in the large limit, through a perturbative analy- 
sis. Therefore it is reasonable to assume that the ground 
state belongs to this sector. 



Ujk — 1 for all links is the obvious choice among the 
configurations of link variables that gives Wp = +l,Vp. 
Of course, any configuration related to this one by a 
gauge transformation will also satisfy the above condi- 
tion on Wp's. With this choice, the Hamiltonian in its 
explicit form becomes 



R 



ci(R)c4(R — as) + JyCi(R)c4(R + ai — as) + JxC2{'R)c3{'R) 

+ JyC3(R)c2(R + aa) + J,ci(R)c2(R) + J3C3(R)c4(R) 

I 



(20) 



where r and a^ are given in Eqs. ([2|) and ([3]). To diago- 
nalizc the Hamiltonian, wc next do a Fourier transform. 



2^ 2^ y_ 27r' 
with /X = 1, 2, 3, 4 and where 

k fcibi + fc2b2 -f fcsbs. 



c=Ma(k), (21) 



and 



b, = ^^, b2^(^^, bs^l 

4 4 2 



Using the property, Cp(— k) = cjj(k), the Hamiltonian 
becomes. 
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H 



dki r dh r dk 



27r 



27r 



2tt 



■4i4(k)c4(k) + 4,4(k)c2(k) + J,ct(k)c2(k) 



(22) 



where Ski = Jx + e '''^^Jy, for i = 1,2. Furthermore, we define Ak — (|<5fciP + + 2J2), and 0k such that 

e~'^^^5kiSk2 = I'^fci I |<5a:2 |e*'^''- The above Hamiltonian is easily diagonahzcd and we obtain the spectrum to be 



i?(k) = ± 



2V2 



{{Jl-\5kA\5u.A) 



' + 2j2(l-cos</)k)|4J|4j} 



(23) 




FIG. 3: In the k space, the contour on which the gap vanishes 
is projected on fca = plane. Jx = Jy = and Jz is varied 
between zero and two. Corresponding values of Jz are shown 
next to the contours. 



In the ground state, all the negative-energy states are 
filled. The system is gapless if solution exists for = 0. 
Since the two terms inside the curly brackets in Eq. (|23p 
are positive definite, for i? = 0, both of them have to 
vanish, i.e., 



J? 



cos (Pk 



+ jy +2 COS 



Jx + Jy +2cOSk2JxJv 



(24) 
(25) 



The values of fci and k2 for which the gap vanishes are 
determined by Eq. (124^ : ^3 is then given by Eq. (pS)) . 
Solutions of Eq. ([24)1 exist only when Jz 1^ Jx + Jy, 



Jx ^ Jy ~t- Jz 



and Jy < Jz 



Jx', these conditions are 



same as that for the 2D Kitaev model. Figure [3] shows 
the plot of contours satisfying Eq. (|24p projected on to 
k^ = plane, where we have set Jx ~ Jy = 1 and varied 
Jz from zero to two. The contour shrinks to the point 



(0,0) as Jz approaches Jx 
opens up. 



Jy = 2, i.e., when the gap 




FIG. 4: The phase diagram: it shows the plane defined 
by Jx + Jy + Jz = Point A corresponds to Jx = 1 and 
Jy = Jz ~ 0, B corresponds to Jy = 1 and Jz = Jx = 0, and 
C corresponds to = 1 and Jx = Jy = 0. The shaded region 
is the gapless phase. 



Figure |4] depicts the phase diagram on a section of 
the parameter phase. Interestingly, the phase diagram is 
symmetric in the three coupling constants though they do 
not appear symmetrically in the Hamiltonian (the z-links 
cannot be transformed to x or y links by any symmetry 
transformation of the lattice) . Nor is the spectrum, given 
by Eq. (j23p . symmetric in the coupling constants. 



We end the discussion on the ground state by noting 
that spin correlations can be calculated in exactly the 
same way as in 2D ^^ < ^^ The only non-zero correlations 
are those constructed out of the interaction terms that 
appear in the Hamiltonian. 



IV. 



FLUX EXCITATIONS 



The Kitaev model in 2D has anyonic excitations, i.e., 
there are particle-like excitations which obey nontrivial 
statistics. Anyons are very specific to 2D and cannot 
exist in higher dimensions; the fundamental reason being 
that in 13 > 2 there are no nontrivial paths which take 
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a particle around another. However, in 3D, excitations 
localized on loops may obey nontrivial statistics. In this 
section, we show that there are excitations in our model 
which are localized on loops. However, we do not address 
the issue of their statistics in this paper. 

We have seen that in the ground state Wp — 1 for 
all plaquettes p. It then follows that the excitations are 
of two types: (1) Flux configurations which violate the 
condition Wp = 1, i.e., Wp = —1 for some of the loops. 
It can be looked upon as creating a tt flux over those 
loops. (2) Fermionic excitations in the background of 
static configurations of Wp. 



We will next show that the excitations of the first type 
have the structure of loops. Earlier we mentioned that 
not all Wp are independent; now we will find the con- 
straints among them. There are four types of plaque- 
ttes which are not related by translation. Let a, h, c, and 
d be the labels for the different types [see Figs. [DJa)- 
EUd)]. Consider a part of the lattice, shown in Fig. [5]Je), 
which consists of four adjacent plaquettes — each a dif- 
ferent type. Here the links are labeled 1-20. Let the 
corresponding loop operators be Wa,Wb, Wc, and Wd, re- 
spectively. In terms of the spin variables these operators 
are 

T-r/ „z „y „^ „x 

Wa = CTiiCri2CTi3Cri4(Ti7Cri8Cr8Cr9CriotTi , 

T-f/ „x X z „y „X „X X 

Wb = CTi4Cri5CrigCr4Cr5Crg(T7Cr8(Tig(Ti7, 

W, = alalalalalalayalalal,, 

Wd = CTnCrj2CTf3Crf4Crf5Cr^6'^|cr|cr|crf. (26) 

Using the relations crjcrjcrl = i and cr"^ — 1, 

WaWbWcWd = 1. (27) 

The above constraint has a graphical interpretation. 
Note that a loop operator can also be written as a prod- 
uct of the interaction terms on the links contained in 
the loop. In Fig. [5l the object (e), which is obtained 
by putting together the four loops (a)-(d), represents the 
left-hand side of Eq. (P?]) . Evidently, each link in (e) 
is shared by two of the Wp's. Since the square of any 
interaction term is one, Eq. (|27p immediately follows. 




FIG. 5: [(a)-(d)] The four types of plaquettes. (e) Part of 
the lattice involving four such adjacent plaquettes; the cor- 
responding operators give rise to a constraint. The ellipses, 
labeled a-d, respectively, represent each of the loops and form 
a "tetrahedron". 



For open boundary conditions, relations such as Eq. ([TT)) 
exhaust all the constraints. To find the configurations 
of {Wp} which are consistent with the constraints, it is 
instructive to consider a lattice obtained by representing 
each plaquette by a single site. A plaquette has a step 
like structure, consisting of two rectangles perpendicular 
to the x-y plane connected by another rectangle on the 
x-y plane. Each loop can be uniquely represented by a 
point at the center of the rectangle on the x-y plane [in 
Fig. [5Jc)such points are marked by ellipses]. Let £ be 



the new lattice thus obtained; topologically, it is the py- 
rochlore lattice, an arrangement of corner-sharing tctra- 
hedra (see Fig. [5]). (The edges of the geometrical object 
formed by the centers of four adjacent loops are not of 
equal length, and hence do not form an exact tetrahe- 
dron but a stretched one. However, the connectivity of 
£ is same as that of the pyrochlore lattice, and we will 
continue to refer to the basic objects as tetrahedron.) 

In this description, the four plaquettes that give rise to 
the constraint in Eq. (P7)) are the four sites of a tetrahe- 
dron, and each tetrahedron corresponds to an indepen- 



Z-llflK 



FIG. 6: The lattice £ formed by the plaquettes — the py- 
rochlore lattice. Loops such as the dashed one, which goes 
through six sites, are the shortest (apart from the basic tri- 
angles) . 



dent constraint. Therefore, any configuration satisfying 
all the constraints will have an even number of Wp taking 
value of —1 in each tetrahedron, where p is now the site 
index in £. Now it is clear how to obtain such config- 
urations: draw a loop C which does not cross itself and 
which lies entirely within the tetrahedra. and let 



Wp = 



-1, ifpeC, 

1, otherwise. 



Any closed self-avoiding loop contains an even number of 
sites (0, 2, or 4) belonging to any particular tetrahedron; 
hence all the constraints are satisfied. In other words, the 
flux excitations have the structure of loops in the lattice 
£,. 



EFFECTIVE HAMILTONIAN IN THE 
LARGE ,h LIMIT 



To study the excitations in the gapped phase, one can 
follow Kitaev'si^ approach to the 2D model — a perturba- 
tive analysis in the limit Jz 3> Jx, Jy The unperturbed 
Hamiltonian Hq and the perturbation Hi are 



Ho = -J. ''^''h (28) 

m = -J. ^^J--^^- E ^f'^l (29) 

<ij>x <ij>y 

H = Ho+Hi. (30) 




FIG. 7: Representation of the flux operators in Heft associated 
with the four types of plaquettes. 



The ground state of Hq is highly degenerate: the two 
spins on each of the z-links can be either | "ft) or | jj, 
). This degeneracy is lifted only at the sixth order of 
perturbation theory. We now calculate the effective low- 
energy Hamiltonian to this order. 

The groundstate subspace of Hq can be thought of as 
the Hilbcrt space of effective spin i's located at each z 
link. Let t^,t^ and be Pauli matrices acting on the 
effective spin at the z link denoted by m, such that, 



TT)m 
-I il) 



(31) 



A plaquette consists of four z-links and six x or y- 
links. Consider a plaquette by also including the two 
z-links that are directed out of the loop (see Fig. [7]). 
Let us label the z-links in such an object with an index 
1 = 1, 2,. ..,6. The labeling scheme for the four different 
types of loops is shown in the figure. We associate with 
each plaquette p a flux operator Bp defined as follows. 



ae 



'pi 'p2 'p3 'pi''pi'''pG ' 



(32) 



where, for i — 2,3,5,6, = a; if the two links in the 
loop attached to the «th z link are the same; and Ui = y 
otherwise. The flux operator associated with the various 
types of loops are also shown in Fig. [71 Then, to sixth 
order — which is the leading order — the effective Hamil- 
tonian, ignoring the constant terms arising at the lower 
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orders, is 



256J5 



(33) 



where the unprimed sum is over the loops of types la and 
2a in Fig. [3 which contain four x hnks and two y hnks, 
while the primed sum is over loops of types 16 and 2b, 
which contain four y links and two x links. Furthermore, 



[Bp, Bp' 



0, 



Bl = 1, yp,p'. 



(34) 
(35) 



Thus, the operators Bp can be simultaneously diagonal- 
ized for all values of p and their eigenstates will also be 
the eigenstates of i?off- From Eq. ([55]) it follows that the 
eigenvalues of Bp are ±1. Therefore, the ground state 
\GS) cs will be such that 



Bp\GS),s^\GS),s, Vp. 



(36) 



It is straightforward to see that Bp is the projection 
of Wp — the conserved quantities of the spin model de- 
fined in Eq. ([5]) — on the ground-state subspace of Hq. 
Therefore, in the large Jz limit, Eq. ([36]) is a confirma- 
tion of our numerically verified assumption earlier that 
the ground state of H is vortex free. 

The excitations arc the states which violate Eq. ([36]) . 
i.e.. Bp = — 1 forp in some subset of plaquettes. However, 
not all Bp are independent — there exist constraints sim- 
ilar to Eq. ((27|) — therefore the plaquettes with Bp = —1 
cannot be chosen arbitrarily. To understand the proper- 
tics of the excitations, such as their statistics, a detailed 
analysis is required. This is currently being done. 



VI. SUMMARY AND DISCUSSION 



fermionic spectrum by mapping the spin model to one 
of free fermions in the background of a static Z2 gauge 
field; the system has a gapped phase and a gaplcss 
one. Quite interestingly, the gap vanishes on a con- 
tour in the k space; this could be related to some ac- 
cidental degeneracies — unrelated to any symmetry of the 
Hamiltonian — of the ground state in the classical limit. 
The two-dimensional Kitacv model has been shown to 
have such a degeneracy in the classical limit, which grows 
exponentially with the system size<^ This result can be 
readily generalized to our model. It will be interesting to 
see how quantum fluctuations lifts this degeneracy. 

We have further shown that the excitations of the 
gauge field, due to some local constraints, have the topol- 
ogy of loops. As a first step towards understanding 
the nature of the excitations in the gapped phase, we 
have derived an effective Z2 gauge theory in the limit 
Jz ^ Jx,Jy The ground state of the effective Hamilto- 
nian thus obtained is trivially solved and it is consistent 
with the assumption we made in the calculation of the 
fermionic spectrum that the ground state is vortex free. 
Further study of the effective theory should tell us more 
about the properties of the excitations, such as whether 
they obey fractional statistics. 

Yao and Kivelson2£ have introduced a two-dimensional 
Kitacv model in which one type of plaquettes is a tri- 
angle; the existence of loops with odd number of links 
results in the spontaneous breaking of time-reversal sym- 
metry. We mention that in our model also the time rever- 
sal symmetry can be similarly broken by replacing each 
vertex in the lattice with a triangle. 

Note added : Toward the completion of this paper, we 
came across the work by Si and Yu in Ref. [H which dis- 
cusses a variety of exactly solvable Kitacv models in three 
dimensions, none of which arc identical to the model in- 
troduced here. 



We have constructed and solved a three-dimensional 
spin-i model which is a generalization of the Kitacv 
model on a honeycomb lattice. Based on the methods 
used by Kitaevji^ we calculated the exact low-energy 



Acknowledgments 

We thank G. Baskaran and R. Shankar for very useful 
discussions. 



* Electronic address: |saptarshi@imsc.res.in| 
t Electronic address: nsurendr@ictp.it 

^ X. G. Wen, Quantum Field Theory of Many-body Systems 
(Oxford University Press, New York 2004). 

2 R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

3 X. G. Wen and Q. Niu, Phys. Rev. B 41, 9377 (1990). 

* G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 

^ G. Baskaran, Z. Zou, and P. W. Anderson, Solid State 
Commun. 63, 973 (1987). 

X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 39, 11413 
(1989). 

N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 (1991). 
^ X. G. Wen, Phys. Rev. B 65, 165113 (2002). 



^ D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61, 
2376 (1988). 

^" N. Read and B. Chakraborty, Phys. Rev. B 40, 7133 
(1989). 

" R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2001). 

M. Freedman, Commun. Math. Phys. 234, 129 (2003). 
P. Shor, Proceedings of the 35th Annual Symposium on 
Fundamentals of Computer Science. (IEEE Press, Los 
Almitos, CA, 1994), pp. 124-134. 
A. Yu. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 
15 A. Yu. Kitaev, Ann. Phys. (N.Y.) 321, 2 (2006). 

C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 



9 



17 


S. D. Sarma, Rev. Mod. Phys. 80, 1083 (2008). 


26 


180404(R) (2007). 


L. M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 


H. Yao and S. A. Kivelson, Phys. Rev. Lett. 99, 247203 




91, 090402 (2003). 


27 


(2007). 


18 


A. Micheli, G. K. Brennen, and P. ZoUer, Nature Phys. 2, 


M. A. Levin and X. G. Wen, Phys. Rev. B 71, 045110 


19 


341 (2006). 


28 


(2005). 


G. Baskaran, S. Mandal, and R. Shankar, Phys. Rev. Lett. 


A. Hamma, P. .Zanardi, and X. G. Wen, Phys. Rev. B 72, 


20 


98, 247201 (2007). 


29 


035307 (2005). 


H. D. Chen and Z. Nussinov, J. Phys. A: Math. Theor. 41, 


H. Bombin, and M. A. Martin-Delgado, Phys. Rev. B 75, 


21 


075001 (2008). 


30 


075103 (2007). 


H. D. Chen and J. Hu, Phys. Rev. B 76, 193101 (2007). 


H. Bombin, and M. A. Martin-Delgado, Phys. Rev. Lett. 


22 


K. P. Schmidt, S. Dusuel, and J. Vidal, Phys. Rev. Lett. 




98, 160502 (2007). 




100, 057208 (2008). 


31 


E. H. Lieb, Phys. Rev. Lett. 73, 2158 (1994). 


23 


S. Dusuel, K. P. Schmidt, and J. Vidal, Phys. Rev. Lett. 


32 


G. Baskaran, D. Sen, and R. Shankar, Phys. Rev. B 78, 


24 


100, 177204 (2008). 


33 


115116 (2008). 


K. Sengupta, D. Sen, and S. Mondal, Phys. Rev. Lett. 100, 


T. Si and Y. Yu, Nucl. Phys. B 803, 428 (2008). 


25 


077204 (2008). 






S. Yang, D. L. Zhou, and C. P. Sun, Phys. Rev. B 76, 







